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ABSTRACT 

Using new and archival radio data, we have measured the proper motion of the black 
hole X-ray binary V404 Cyg to be 9.2 ± 0.3masyr~^. Combined with the systemic 
radial velocity from the literature, we derive the full three-dimensional heliocentric 
space velocity of the system, which we use to calculate a peculiar velocity in the range 
47-102 km s~^, with a best fitting value of 64kms~^. We consider possible explanations 
for the observed peculiar velocity, and find that the black hole cannot have formed 
via direct collapse. A natal supernova is required, in which either significant mass 
(~ IIMq) was lost, giving rise to a symmetric Blaauw kick of up to ^ 65kms^^, 
or, more probably, asymmetries in the supernova led to an additional kick out of 
the orbital plane of the binary system. In the case of a purely symmetric kick, the 
black hole must have been formed with a mass ~ 9-^©, since when it has accreted 
0.5-1.5 Mq from its companion. 

Key words: X-rays: binaries - astrometry - radio continuum: stars - stars: individual 
(V404 Cyg) - stars: supernovae: general - stars: kinematics 



and distance to the system, the full three-dimensional space 
velocity of the system can be derived. Along with system 
parameters such as the component masses, orbital period, 
donor temperature and luminosity, these parameters may 
be used to reconstruct the full evolutionary history of the 
binary system back to the time of compact object forma- 
tion, as done for the systems GRD_jJ_655-40_^^niemr^er^l] 
l2005h and XTE J 1118-^480 (|Fragos et al.ll2007l ). 

By studying the distribution of black hole X-ray binary 
velocities with compact object masses, we can derive con- 
straints on theoreti cal models of black hole formation (e.g. 
iFrver fc Kalogcr a 20011 ) . The two most common theoretical 
scenarios for creating a black hole involve either a massive 
star collapsing directly into a black hole with very little or 
no mass ejection, or delayed formation in a supernova, as 
fallback onto the neutron star of material ejected during 
the explosion creates the black hole. In the latter case, con- 



1 INTRODUCTION 

The proper motions of X-ray binary systems can be used 
to derive important information on the birthplaces and for- 
mation mechanisms of their compact objects. Since typical 
proper motions are of order a few milliarcseconds per year, 
high-resolution observations and long time baselines are re- 
quired to measure the transverse motions of such systems 
across the sky. Proper motions have only bee n measured for 
a handful of X-ray binary systems to date ([M irabel et al.l 
2001, 2 0(]^:lRib6 et al.ll2q02l : [Mirabel fc Rodrigi ics 2003a-h: 
bhawan^et al.ll2006l . |2007^ ■ and only one X-ray binary, Sco 



X-1, has a meas ured proper motion, parallactic distance, and 
radia l velocity (|Bradshaw et al.l Il999l : ICowlev fc CramptonI 
1 19751 ). With the position, proper motion, radial velocity. 
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straints on the magnitude and symmetry of any natal kick 
can also be derived. When the primary star reaches the end 
of its life and explodes as a supernova, the centre of mass of 
the ejected material continues to move with the velocity the 
progenitor had immediately before the explosion. The cen- 
tre of mass of the binary system then r ecoils in the opposite 
direction. Such a kick (|Blaauwl [l96ll ) is thus constrained 
to lie in the orbital plane. In addition, in the presence of 
asymmetries in the supernova explosion, a further, asym- 
metric, kick, which need not be in the orbital plane, may be 
impa r ted to the binary (see, e.g., Brandt fc Podsiadlowskil 



1 19951 : iPortegies Zwart fc YungelsonI Il998l : iLai et all I2OOII . 
for more detailed over views). While pulsar proper motions 
l|Lvne fc Lorimeilll994l ) and the high eccentricities of Be/X- 
ray binaries provide strong evidence for asymmetric kicks 
in neutron star formation dvan den Heuvel fc van ParadiisI 



I1997I : Ivan den Heuvel et al.H2000l ). less attention has been 
paid to black hole binary systems. In two of the black hole 
X-ray binaries with some of the best observational con- 
straints on the system parameters, XTE J 1118-1-480 and 
GRO J1655-4 0, there is evidence for an asymmetric su - 
pernova kick (jGualandris et al.ll2005l : IWillems et al.ll2005h . 
There is also some evidence suggesting an asymmetric kick 
in GRS 1915-1-105 (|Dhawan et al.l l2007h . although in that 
case the uncertainty in the distance to the system pre- 
cludes the derivation of strong constraints on the kick ve- 
locity. However, the black hole in C ygnus X-1 is inferred 
to hav e formed via direct collapse (jMirabel fc RodriguesI 
l2003al ). More black hole sources need to be studied to mea- 
sure black hole kicks and investigate consequent formation 
mechanisms in order to constrain theoretical models of black 
hole formation. 



1.1 V404 Cyg 

V404 Cyg is a dynamically-confirmed black hole X-ray 
binary system, with a m ass function of 6.08 ± 0.06Mq 
lluasares fc Charledll994h . The system comprises a black 
hole accretor of mass 1212-^0 (|Shahbaz et al.l [19941 ) in a 
6.5-d orbit with & Q.ltnlMo) KO IV stripped giant star 
(|Casares et al.lll993l : lKing||l993l '). The orbit is highly circular, 
with an eccentricity of e < 3 x 10"*, in agreement with the 
fact that tidal forces just before the onset of mass transfer 
should have recircularised the orbit follo wing the supernova. 
Its lo w radial velocity (— 0.4±2.2 kms~^: [Casares fc Charlea 
1 19941) has been taken as evidence for a small natal kick 
l|Brandt et al.. 1995; Nelemans et al.lll999l ). 

In this paper we present measurements of the proper 
motion of V404 Cyg, and go on to derive its full three- 
dimensional space velocity, infer its Galactocentric orbit, 
and discuss the implications for the formation of the black 
hole in the light of the inferred natal kick from the super- 
nova. 



2 OBSERVATIONS AND DATA REDUCTION 

In order to investigate the proper motion of the source, 
we interrogated the Very Large Array (VLA) archives for 
high-resolution observations of V404 Cyg. We selected only 
A-configuration observations at frequencies of 8.4 GHz and 
higher, to obtain the highest possible astrometric accuracy. 



We further restricted the dataset to observations where 
the phase calibrator was J 2025-f 3343, the phase reference 
source u sed in the High Sensitivi ty Array (HSA) observa- 
tions of iMiller- Jones et al] (|2008l ). Since all positions are 
measured relative to the phase reference source, observations 
using a different secondary calibrator would potentially have 
been subject to a systematic positional offset. 



The VLA data were reduced using s tandard proce- 
dures within the 31Dec08 version of MPS (|Greisenll2003l ). 
A script was written in ParselTongue, the Python in- 
terface to AlPS, to automate the bulk of the calibration. 
We corrected the derived positions for shifts in the as- 
sumed calibrator position, using a reference position of 
20''25™10:8421050 33°43'00'.'214430 (J 2000) for the calibra- 
tor source. All co-ordinates were precessed to J 2000 val- 
ues using the AlPS task UVFIX. Using data from the USNO 
(http:/ /maia. usno.n avy.mil/ser7/finals2000 A. all J , we cor- 
rected for offsets in (UTl-UTC), the difference between Uni- 
versal Time (UTl) as set by the rotation of the Earth and 
measured by very long baseline interferometry (VLBI), and 
co-ordinated Universal Time (UTC), the atomic time (TAI) 
adjusted with leap seconds. Where measured offsets were 
available, we also corrected for shifts in antenna positions 
using the AlPS task vlant. Source positions were measured 
by fitting an elliptical Gaussian to the source in the image 
plane, using the deconvolved, phase-referenced image prior 
to any self-calibration. The source was not resolved in any 
of the images. We added an extra positional uncertainty of 
10 mas to the measured VLA positions, to account for sys- 
tematic uncertainties in the astrometry. The list of observa- 
tions and derived source positions is given in Table [T] 



The dataset was enhanced by the use of two high- 
resolution VLBI measurements of the source position . We 
used the 8.4-GHz position of Miller- Jones et ahl (|2008t ). and 
also obtained a second measurement using global VLBI at 
22 GHz, under proposal code GM064. Eight European sta- 
tions (Cambridge, Effelsberg, Jodrell Bank Mkll, Medicina, 
Metsahovi, Noto, Onsala, and Robledo), all ten Very Long 
Baseline Array (VLBA) stations, the phased VLA and the 
Green Bank Telescope (GBT) participated in the experi- 
ment. Data were taken from 21:30:00 UT on 2008 May 31 
until 16:00:00 UT on 2008 June 1, with the European stations 
being on source for the first 12 h of the run (Robledo from 
02:20:00 until 08:45:00 on 2008 June 1) and the North Amer- 
ican stations for the second 12 h (Mauna Kea from 07:30:00, 
when the source rose in Hawaii). The overlap time, when 
both sets of stations were on source, was 5.5 h. The phase 
reference and fringe finder source was J 2025-1-3343. We ob- 
served with a total bit rate of 512 Mb s~^ , with a bandwidth 
of 64 MHz per polarization. We observed in 150-s cycles, 
spending 1.5 min on the target and Imin on the calibra- 
tor in each cycle. The VLA was phased up at the start of 
each calibrator scan, and we made referenced pointing ob- 
servations with the larger dishes (the VLA, GBT, Effelsberg 
and Robledo) every 1-2 h. Data were reduced using standard 
procedures within AlPS. We detected the source at a signifi- 
cance level of 4.7cr. 
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Code 


Date 


MJD 




Config. 


Frequency 
(GHz) 




R.A 


Error 
(sec) 


Dec 




Error 

( 3,rCSGC J 


AH348 


1990 Fob 16 


47938.71 


0.18 


A 


14.94 


20''24 


"03^^82854 


0.00085 


33°52 


02^ 


0348 


0.0103 


AH385 


1990 Mar 08 


47958.70 


0.12 


A 


14.94 


20** 24 


"03? 82896 


0.00083 


33° 52 


02^ 


0379 


0.0102 


AH390 


1990 Mar 25 


47975.55 


0.09 


A 


8.44 


20** 24 


"03? 82865 


0.00081 


33°52 


02^ 


0439 


0.0101 


AH424 


1991 Sop 25 


48525.03 


0.05 


AB 


14.94 


20''24 


"03^82733 


0.00276 


33° 52 


01^ 


9563 


0.0235 


AH424 


1991 Sop 25 


48525.04 


0.05 


AB 


8.44 


20'>24 


"03^82321 


0.00148 


33° 52 


02^ 


0143 


0.0136 


AH424 


1992 Oct 20 


48916.03 


0.07 


A 


14.94 


20''24 


"03? 82752 


0.00096 


33° 52 


02^ 


0208 


0.0122 


AH424 


1992 Oct 20 


48916.04 


0.07 


A 


8.44 


20''24 


"03=82673 


0.00084 


33° 52 


02^ 


0196 


0.0105 


AH390 


1993 Jan 28 


49015.68 


0.04 


AB 


14.94 


20''24 


"03? 82779 


0.00222 


33° 52 


01^ 


9964 


0.0173 


AH390 


1993 Jan 28 


49015.68 


0.02 


AB 


8.44 


20''24 


"03? 82427 


0.00219 


33° 52 


02^ 


0317 


0.0158 


AH641 


1998 May 04 


50937.60 


0.04 


A 


8.46 


20'' 24 


"03? 82493 


0.00086 


33° 52 


01^ 


9820 


0.0109 


AH669 


1999 Jul 04 


51363.41 


0.01 


A 


8.46 


20''24 


"03? 82479 


0.00144 


33° 52 


01^ 


9721 


0.0156 


AH669 


1999 Jul 13 


51372.37 


0.01 


A 


8.46 


20'>24 


"03?82431 


0.00115 


33° 52 


01^ 


9984 


0.0143 


AH669 


1999 Jul 26 


51385.30 


0.01 


A 


8.46 


20'>24 


"03? 82508 


0.00260 


33° 52 


01^ 


9470 


0.0343 


AH669 


1999 Aug 22 


51412.18 


0.01 


A 


8.46 


20'' 24 


"03=82443 


0.00163 


33° 52 


01^ 


9429 


0.0176 


AH669 


1999 Sep 01 


51422.21 


0.01 


A 


8.46 


20'' 24 


"03? 82417 


0.00152 


33° 52 


02^ 


0063 


0.0207 


AH669 


1999 Sop 04 


51426.02 


0.01 


A 


8.46 


20'' 24 


"03?82397 


0.00194 


33° 52 


01^ 


9464 


0.0188 


AH669 


2000 Oct 20 


51837.33 


0.01 


A 


8.46 


20'>24 


"03=82718 


0.00595 


33° 52 


01^ 


9676 


0.0405 


AR476 


2002 Fob 03 


52308.69 


0.01 


A 


8.46 


20'' 24 


"03? 82325 


0.00088 


33° 52 


01^ 


9481 


0.0108 


AR476 


2002 Mar 01 


52334.54 


0.01 


A 


8.46 


20'>24 


"03=82586 


0.00257 


33° 52 


01^ 


9457 


0.0224 


AH823 


2003 Jul 29 


52849.30 


0.29 


A 


8.46 


20'' 24 


"03? 82440 


0.00101 


33° 52 


01^ 


9317 


0.0120 


BG168 


2007 Dec 02 


54436.84 


0.09 


HSA 


8.42 


20'' 24 


"03=82129 


0.000010 


33° 52 


01^ 


8993 


0.0003 


GM064 


2008 May 31 


54618.42 


0.25 


Global 


22.22 


20'' 24" 


^03? 82 1082 


0.000009 


33° 52 


01^ 


8957 


0.0001 



Table 1. Summary of the VLA and VLBI observations. The errors on the MJD values are taken as half the length of the observation. The 
quoted positional errors for the VLA observations are the sum in quadrature of the statistical errors and a systematic uncertainty of 10 mas. 
All co-ordinates are for epoch J2000. The coordinate system is based on the assumed position of the calibrator source J 2025+3343, taken 
to be (J 2000) 20''25™10!8421050(224) 33°43'00'.'21443(42). The HSA and global VLBI observations have been corrected for parallax, 
assuming a source distance of 4kpc, and accounting for the uncertainty in the parallax when calculating the positional error bars. 



3 RESULTS 

Fig. [T] shows the measured Right Ascension and Declination 
as a function of time, over the ~ 20 years since the 1989 
outburst of V404 Cyg. The best fitting proper motions, in 
Right Ascension and Declination respectively, are 

cos 5 = -4.99 ±0.19 mas y"^ (1) 
/ii- = -7.76 ± 0.21 mas y"\ (2) 

Thus the total proper motion is jj. — 9.2 ± O.Smasy^^. All 
uncertainties are 68 per cent confidence limits, and unless 
otherwise noted we will henceforth quote la error bars on 
all measurements. 

The fit gives a reference position (prior to correcting for 
the effects of the unknown parallax) of 20''24'"03';82177(2) 
33°52'01'.' 9088(3) (J 2000) on MJD 54000.0, from which we 
can use the proper motion to determine the predicted posi- 
tion at any future time. 



3.1 Converting to Galactic Space Velocity 
co-ordinates 

With the systemic radial velocity of —0.4 ± 2.2kms~^ 
measured from optical Hq studies (|Casares fc Charles! 
1 19941 ). we can calculate, for a given source distance, 
the full three-dimensional space velocity of the system. 
The best constraint on the source distance is 4.0^?'o kpc 

I 1 I 1 —l.Z r- 

(IJonker fc Nclcmans 2004). Using the transformations of 
IJohnson fc Sodcrblom (198^), and the standard solar mo- 
tion of {Uq = 10. ± 0.36, Fp) = 5 25 ± 0.62, We = 
7.17 ±0.38) kms"^ jDehnen fc Binnevl[l99i ). we can com- 
pute the heliocentric Galactic space velocity components U, 
V and W (defined as U positive towards the Galactic Cen- 
tre, V positive towards / = 90°, and W positive towards the 
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Figure 1. Top panel: Measured Right Ascension as a function of 
time, from 1990 to 2007. Bottom panel: Declination as a function 
of time. The dotted lines are the best fitting proper motions, 
= (-4.99 ±0.19) X 10"''secy-i in R.A. and fig = (-7.76 ± 
0.21) masy~l in Dec. The black dots are the VLA 8.4 GHz points 
and the light grey triangles are the 15-GHz VLA measurements. 
The dark grey square s labelled 'H' and 'G' are from the 8.4-GHz 
HSA observations of iMiller- Jones et al.l |20oi) and the 22-GHz 
global VLBI observations reported in this paper, respectively. 



North Galactic Pole). The derived system parameters are 
given in Table O 

For a given distance, the expected values of U and V 
can be calculated , assuming t he source participates in the 
Galactic rotation. iReid fc Brunthal er (.2004,) determined the 
angular rotation rate of the LSR at the Sun, Qq/Rq — 
29.45 ± 0.15kms ~^ kpc~^. A ssuming a Galactocentric dis- 
tance of 8.0 kpc (jReidlll993h . this implies a circular veloc- 
ity of 236kms~^. For circular rotation about the Galactic 
Centre, the W component of the velocity is expected to be 
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Parameter 



Value 



Galactic longitude I 

Galactic latitude b 

Distance d (kpc) 

Systemic velocity 7 (kms^-"^) 

Proper motion /Iq cos 5 (masy^^) 

Proper motion fi^ (masy~^) 

U (kms-l) 

V (kms-l) 

W (kms-l) 

U- <U > (kms-i) 

V- <V > (kms-i) 

W- <W > (kms-i) 

Upec (kms~^) 



73.12° 
-2.09° 
4.0 



+2.0 
1.2 
0.4 ±2.2 



-4.99 ±0.19 

-7.76 ±0.21 

177.1 ±3.8i«g-i 

-46.1 ± 2Atfd 

0.2 ±3.7135 

62.0lf,l ■ 
1+6.5 
-0.0 



-16.1"' 



64.1±3.7tigg 



Table 2. Measured and derived parameters. U, V and W are 
the Galactic space velocity components in the direction of the 
Galactic Centre, / = 90° and b = 90° respectively. The first set 
of error bars accounts for uncertainties in the measured space 
velocities only, and the second takes into account the distance 
uncertainty. U- <U >, V- <V > and W- < > are the dis- 
crepancies from the velocities that would be expected for Galactic 
rotation. Summing these discrepancies in quadrature gives the pe- 
culiar velocity, Vpcc- The major source of error in these values is 
the distance uncertainty. 



zero. As done by iDhawan et all (|2007l ) for GRS 1915+105, 
we can transform the measured values of U and V into radial 
and circular velocities about the Galactic Centre, expected 
to be Drad = Okms^^ and Ucirc = 236kms~^ respectively. 
Fig. [2] shows the derived radial, circular and W velocities 
and the peculiar velocity of V404 Cyg. The peculiar veloc- 
ity is defined to be the difference between the measured 3- 
dimensional space velocity and that expected for a source 
participating in the Galactic rotation, 

«pcc = («?ad + {Vci.c - 236)' + W^) ^'^ (3) 

F or the range of distances found bv lJonker fc NelemansI 
(|2004l ). 2.8-6.0 kpc, it is clear that the peculiar velocity is 
non-zero. We derive a value of Upec = 64.1 ± 3.7ljg g kms~^, 
where the first error bar accounts for statistical error in the 
space velocities, and the second for the distance uncertainty. 
The predominant component of the peculiar velocity is ra- 
dial, with a circular velocity slightly faster than expected, 
and a velocity out of the Galactic plane consistent with zero. 



4 THE ORBITAL TRAJECTORY IN THE 
GALACTIC POTENTIAL 

From the known source position and the measured spatial 
velocity components, we can integrate backwards in time 
to compute the orbital trajectory of the system in the po- 
tential of the Galaxy. Usin g a fifth-order Runge-Kutta al- 
gorithm (|Press et al.lll993 ) to perform the integration, we 
compared the predictions of several different models for the 
Galactic potential (C arlberg fc Innanen 1987, u s ing the re- 
vised parameters of Kuiiken fc Gilmord 19891: Paczvhski 



I990I : I Johnston et al.|[r995l; IWolfire et al.H 19951 : fFlvnn et al 
19961 : de Oliveira et al. 120021 ). all using some combination 
of one or more disc, spherical bulge and halo components 
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Figure 2. Derived Galactoccntric velocities of V404 Cyg, as a 
function of source distance, (a) shows the radial velocity, Uradi 
(b) shows the circular velocity, iicirci (c) shows the velocity out of 
the Galactic plane, W , and (d) shows the peculiar velocity, Vpcc, 
i.e. the difference from the values expected for a source partici- 
pating in the Galactic rotation. Dotted lines show the expected 
values of 0kms~^ for the radial, W and peculiar velocities, and 
236 km s^^ for the circular velocity. Error bars only account for 
uncertainties in the space velocity components, assuming zero er- 
ror on the distance at each plotted point. 



A representative orb it reconstruction using the model of 
Ijohnston et al.l (|l995h is shown in Fig. [3] While the errors 
in the space velocity components make little difference to 
the computed orbital trajectory for a given model, the un- 
certainty in the distance has much more of an effect. Fig. |4] 
compares the trajectories computed for 2.8, 4.0 and 6.0 kpc, 
the lower, mean and upper bounds to the possible range of 
distances. A larger source distance implies a more ellipti- 
cal orbit, which reaches further from the Galactic Centre at 
apogalacticon. 

Comparing the different models for the Galactic poten- 
tial, we find that the predicted orbital trajectories in the 
Galactic plane begin to diverge significantly after only 25- 
30Myr. In the perpendicular direction, the predictions di- 
verge even faster, within 2-3 Myr. Given the uncertainty in 
the source distance and Galactic potential, it is clearly im- 
possible to integrate back in time over 0.4-0.8 Gyr (Section 
O to locate the birthplace of the system. However, for a 
given source distance, we can average over the ensemble of 
model predictions to derive some generic properties of the 
orbit. The eccentricity e, semi-major axis a, and the dis- 
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Figure 3. The computed orbit of V404 Cyg over the last 1 Gyr, using the Galactic potential of I Johnston et al] lll995l ) and assuming 
a current source distance of 4.0 kpc. In (a), (b) and (c), the greyscale is a logarithmic representation of the Galactic mass density 
(Mq kpc~^), the black line denotes the computed orbit for the best fitting space velocities and positions, and the white trace indicates 
the error due to uncertainty in the space velocity components. For clarity, the effect of the distance uncertainty is not shown (see Fig. |4)l. 
Panels (d) and (e) are a zoomed-in version of the x—z and y—z planes respectively (x, y and z defined as the directions I = 0° , I = 90°, 
b = 90° respectively), with the black line indicating the best-fitting trajectory and the grey trace once again indicating the spread due 
to uncertainty in the measured space velocities only. In all panels, the cross marks the Galactic Centre, and the open circle and triangle 
mark the current positions of the Sun and V404 Cyg respectively. 



tances of the apsides, Tmax and Tmin, for the motion in the 
Galactic plane, the maximum height reached above or be- 
low the plane, Zmax, and the minimum and maximum values 
of the peculiar velocity, Wpcc.min and Wpoc.max, are given in 
Table El 

We can compare thes e values with those of the thick 
and thin disk populations. iBinnev fc Merrifieldl (|l998l ) give 
the velocity dispersions in the radial, azimuthal and vertical 
directions (relative to the Galact ic plane) for both stella r 
populations, based on the data of lEdvardsson et al.l (|l993l ). 
for stars within 80 pc of the Sun. A set of Monte Carlo sim- 
ulations established that the typical planar eccentricities of 
the thin and thick disk populations were 0.12 ± 0.06 and 
0.29 ± 0.15 respectively, while the typical maximum height 
reached above the plane was 0.17 ± 0.14 and 0.56 ± 0.61 kpc 
respectively. Unless V404 Cyg is at the maximum possible 
distance, it is likely to have originated as a member of the 
thin disk population (as also suggested by the age, metallic- 



ity and component masses of the system), and received some 
sort of kick which increased the planar eccentricity and the 
component of velocity out of the Galactic plane. 



5 THE PECULIAR VELOCITY 

The current peculiar velocity of the system is 64.1 ± 
3.7lig g kms~^. However, owing to its orbit in the Galac- 
tic potential, this is not a conserved quantity. For a distance 
of 4 kpc, we find that the peculiar velocity varies between 39 
and 79kms"^ (see Tableland Fig.[5l). We go on to examine 
potential explanations for this peculiar velocity. 

5.1 Symmetric supernova kick 

If it formed with a natal sup ernova, the system can receive 
a Blaauw kick (|Blaauw.il96ll ) , whereby the binary recoils to 



© 2008 RAS, MNRAS OOClTlfTOl 



6 J.C.A. Miller- Jones et al. 





2.8kpc 




— 4.0kpc 




— e.Okpc 




-20 -15 -10 -5 5 
y(kpc) 



10 15 20 



Figure 4. The effect of the distance uncertainty on the com- 
puted Galactocentric orbit of V404 Cyg. Trajectories have been 
computed for the minimum, be s t fitti ng and maximum distances 
derived bv ljonker &: NelemansI ||2004| ) , using the measured s pace 
velocities and the Galactic potential of ljohnston et all |l99^ . For 
each distance, the spread in the trajectories due to uncertainties 
in the space velocity components has been plotted. In all panels, 
the cross marks the Galactic Centre, and the open circle and tri- 
angles mark the current positions of the Sun and V404 Cyg (at 
its different distances) respectively. 



conserve momentum after mass is instantaneously ejected 
from the primary. For the binary to remain bound after the 
supernova, the ejected mass AM must be less than half the 
total mass of the system. A maximum ejected mass trans- 
lates to a maximum recoil ve locity of the binary, fo r which 
an expression was derived bv lNelemans et al] (|l999l l. 



m 



Pro 



d 

Mbh -I- m 



-1/3 



M, 







-5/3 



kms 



(4) 



where Mbh and m are the masses of the black hole and 
the secondary immediately after the supernova, before mass 
transfer has begun, and Pro-circ is the period of the or- 
bit once it has recircularized following the supernova, re- 



Distance 



2.8 kpc 



4.0 kpc 



1.0 kpc 



rmax (kpc) 


9.8 ±0.4 


11.7 ±0.7 


20.2 ± 4.4 


rmin (kpc) 


7.0 ±0.1 


7.2 ± 0.1 


7.9 ±0.1 


2max (kpc) 


0.13 ±0.02 


0.20 ±0.03 


0.47 ±0.14 


a (kpc) 


8.4 ±0.2 


9.4 ±0.4 


14.1 ± 2.2 


e 


0.16 ±0.02 


0.24 ±0.03 


0.42 ±0.07 


^'pcc,miii (kms^^) 


21.9 ±5.1 


38.9 ±6.4 


82.4 ±7.4 


^pcc,max (km S ) 


56.5 ±4.9 


78.8 ±5.7 


134.3 ± 10.3 



Table 3. Derived parameters of the Galactocentric orbit, av- 
eraged over the ensemble of models for the Galactic potential. 
The parameters are the maximum and minimum distances from 
the Galactic Centre measured in the plane (rmax and r^^^), the 
maximum distance reached above or below the plane (2:niax 

), the 

semi-major axis a, the orbital eccentricity e (both calculated in 
the plane), and the minimum and maximum peculiar velocity 
(^pec.min ^^nd Dpcc,max respectively). Uncertainties are the scat- 
ter due to both the errors on the measured velocities and the 
differing models for the Galactic Potential. 



lated to the orbital period immediately after the supernova 



has occurred by Pr 



re — circ 



= P 



post- 



-sn(1 



-post— SNy 



N3/2 



Re- 



circularization occurs before the start of mass transfer from 
the secondary to the black hole, due to the strong tidal forces 
present once the secondary has evolved and expanded suffi- 
ciently to come close to filling its Roche Lobe. 

During its X-ray binary phase, the system undergoes 
mass transfer from the donor star to the black hole, increas- 
ing the black hole mass and orbital period and reducing the 
donor mass. In the case of conservative mass transfer, angu- 
lar momentum conservation and Kepler's Third Law imply 
that these parameters evolve as 



Pcx {M^nmy 



(5) 



While we have no exact constraint on how long ago mass 
transfer began, we ca n use the current orbital period 
and component masses (|Shahbaz et al.lll994l ) together with 
Equation [5] to find the system parameters at any point in 
the past, and the implied maximum ejected mass and recoil 
velocity of a supernova which would have created a system 
with those parameters (from Equation 3)) . This is shown in 
Fig. [6] for three possible current configurations of donor and 
accretor mass. 

As mass is transferred from donor to accretor, the or- 
bital period increases (Equation (S}. Thus for a larger trans- 
ferred mass, Mtrans, (a smaller initial black hole mass and 
a less recent onset of mass transfer) the orbital period at 
the onset of mass transfer is smaller. The pre-supernova or- 
bital period, which is smaller due to the mass loss in the 
explosion (middle line in the top left hand panel of Fig. 
cannot be so short that either the black hole progenitor or 
its companion fills its Roche lobe, setting a minimum pos- 
sible pre-supernova orbital period (lower line in the top left 
hand panel of Fig. (Gjl. For small values of Aftrans, the max- 
imum mass lost in the supernova is equal to half the total 
system mass, and the maximum kick velocity Umax increases 
via Equation |4] because both the post-supernova donor mass 
increases and the post-supernova period decreases with in- 
creasing Mtrans- Where the period constraint becomes im- 
portant, the maximum mass lost in the supernova decreases, 
being set by the shortest allowed pre-supernova orbital pe- 
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Figure 5. Variation in Galactocentric radial velocity i;radi circu- 
lar velocity iJcirci perpendicular velocity Vz, and peculiar velocity 
Vpec, as a function of time while integrati ng the orbit backward s 
over the last Gyr, using the potential of I Johnston et all |l99^. 
Solid, dotted and dashed lines are for source distances of 4.0, 
2.8 and 6.0 kpc respectively. The grey lines show the uncertainty 
arising from the error bars on the measured space velocities, for 
a distance of 4.0 kpc. The peculiar velocity varies as a function of 
time. 




Figure 6. Evolution of the post-supernova binary orbital pe- 
riod (top panels), maximum possible mass ejected in the super- 
nova (middle panels), and maximum possible velocity kick (lower 
panels), as a function of total mass transferred since the super- 
nova from the secondary to the black hole. Left-hand panels are 
all for the case of conservative mass transfer, for current com- 
ponent masses of {M^ii,m) = (12, 0.7) Mq. The top left plot 
shows the post-supernova orbital period (thick solid line) , the pre- 
supernova orbital period (middle line) and the minimum permit- 
ted pre-supcrnova orbital period in which neither the helium star 
progenitor nor the main sequence companion fills its Roche lobe 
(bottom line). The right-hand panels show the situation for con- 
servative mass transfer using different current component masses, 
and one case of non-conservative mass transfer. Thick solid lines 
are for current component masses of {M^ii,m) = (12, 0.7) Mq, 
thin solid lines are for (10, 1) Mq and grey lines are for (14, 0.5) 
Mq, assuming conservative mass transfer in all cases. Dashed 
lines indicate non-conservative mass transfer for the (10, 1) Mq 
case, with ten per cent of the mass transferred being lost from 
the system. 



riod. The decreasing mass loss in the explosion then offsets 
the increasing post- supernova donor mass, and the max- 
imum possible kick velocity decreases as Mtrans increases 
further. 

Equation[5]is valid only for the case of conservative mass 
transfer. However, there are strong indications that this as- 
sumption is not valid for V404 Cyg. At the very least, the 
radio outbursts (e.g. iHan fc Hiellming 1 19921) suggest that 
mater ial is being lost to jet outflows. Podsiadlowski et al.l 
(|2003l ) calculated binary evolution tracks for a WMq black 
hole in orbit with donors of mass 2-17Mq (for there to 
have been any symmetric kick, Fig. [6] shows that the trans- 
ferred mass must be < 2.5Mq, implying that the initial sec- 
ondary must have been less massive than ~ 3.5Mq), some 
of which they found could reproduce the system parame- 



ters of V404 Cyg extremely well. They took into account 
non-conservative mass transfer, whereby mass transferred 
in excess of the Eddington rate is lost from the system. A 
comparison with the equations for conservative mass trans- 
fer suggested that the orbital period increases more slowly 
in the non-conservative case, by a factor of at most 2. As an 
example of how this could affect the maximum Blaauw kick 
velocity, the dashed lines in Fig. [6] are for the current case of 
a IOMq accretor with a IMq donor, with 10 per cent of the 
transferred mass being lost from the system and the orbital 
period increasing a factor 2 more slowly than predicted by 
Equation [5] 

Our calculations show that for V404 Cyg, it is just pos- 
sible to achieve a kick of 64kms~'^ using only symmetric 
mass loss in the supernova (a Blaauw kick). If the explo- 
sion occurred at a point in the Galactocentric orbit where 
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the peculiar velocity was minimized (Fig. [Sjl , this could be 
sufficient to explain the observed peculiar velocity, particu- 
larly if the system is towards the lower end of the possible 
range of distances. However, to achieve such a large kick re- 
quires fine-tuning the parameters, to give a relatively low 
initial black hole mass (closer to Q.2Mq than the best fit- 
ting value of 12Mq) and an ejection of ~ IIM© during 
the supernova. Such a large mass loss, in combination with 
the system parameters (a ~ 9Mq black hole with a low- 
mas s donor), do not make t his scenario very plausible (see, 
e.g.. lFrver fc Kalogerall200ll . for theoretical estimates of the 
amount of mass ejected in a supernova as a function of pro- 
genitor mass; to form a QMq black hole requires a progenitor 
of 25-30 M0, in which case the mass ejected in the super- 
nova is expected to be < 3Mq). We find it unlikely that 
a symmetric kick alone is sufficient to explain the observed 
peculiar velocity. 



5.2 Velocity dispersion in the disc 

The donor star in V40 4 Cyg is a KO subgiant 
( Casares fc Charled 1 19941') of mass 0.7lo:2A/o 
( Shahbaz et all Il994l '). iKind (|l993l ) demonstrated that 
the system is a stripped giant, which evolves on the 
nuclear timescale of the donor star, which is in the range 
0.4-0.8 Gyr. Thus the secondary was initially significantly 
more massive, and has transferred mass to the black hole. 
From Fig. |6l the donor star is unlikely to have transferred 
more than 2.5M0 to the black hole during the mass transfer 
phase, implying an i nitial donor mass of < S .SMp. The 
evolutionary tracks of lPodsiadlowski et al.l (|2003l ) show that 
such a system would take of order 0.7-0.8 Gyr to evolve 
to an orb ital p e riod o f 6.5 d, in agreement with the range 
given bv iKin j (|l993l ). In 0.8 Gyr, the system has made 
3.5-5 orbits in the potential of the Galaxy (depending on 
the model used for the Galactic potential; see Section |4]), 
so could well have received some component of peculiar 
velocity in the Galactic plane owing to non-axisymmetric 
forces such as scattering from the potentials of spiral arms 
or interstellar cloud complexes l|Wielenl Il977l ). However, 
estimat es of the velocity disper s ion of the thin disk popu- 
lation (|Dehnen fc BinnevI Il998l : iMignardI l2000h show that 
for F0-F5 stars (with initial masses comparable to that 
estimated for the donor star in V404 Cyg prior to the onset 
of mass trans f er), it is of order 28kms~^. 

iMignardI (|2000l ) fitted the Sun's peculiar motion and 
the differential velocity field caused by the Galactic rota- 
tion to the parallaxes and proper motions measured by Hip- 
parcos. Examining the fit residuals showed the peculiar ve- 
locities to follow a three-dimensional Gaussian distribution, 
with the predicted scatter in the velocity in Galactic lon- 
gitude being given by (0.7 < >)^^^, where < > is 
the velocity dispersion in the radial direction, determined as 
< >i/2= 22.5±0.3kms"^ for F0-F5 stars. The measured 
peculiar velocity of V404 Cyg in the longitudinal direction 
is — 64.0137 7 km s~^, implying a probability of 7 x 10^"' of 
being caused by Galactic velocity diffusion if the source is 
at 4kpc, and only 1.6 x 10~^ even if the source is at 2.8 kpc. 
Thus this mechanism is unlikely to account for the observed 
and inferred range of peculiar velocities. 



5.3 An asymmetric supernova kick? 

If neither the effects of stellar diffusion nor a symmetric 
kick can explain the observed peculiar velocity, it could in- 
stead be the effect of an asymmetric kick during the super- 
nova. The smaller dispersion in ^-distance (distance above 
or below the Galactic plane) of black hole X-ray binaries 
when compared to neutron star systems had been inter- 
preted as eviden ce for smaller kicks whe n forming black holes 
( White fc van Paradiis .1996.). However. I Jonker fc NelemansI 
( 20041 ). using a larger sample and revised distance estimates 
for the sources, found no evidence for such a discrepancy, 
suggesting that black holes can receive natal kicks com- 
parable with those seen in neutron star systems. Recent 
analysis of the space vel ocities of the black hol e X-ray bi- 
naries XT E J 1118-H480 dCualandris et al.ll2005l ) and GRO 
J 1655-40 (IWiUems et al.ll2005l ) found evidence for asymmet- 
ric kicks in these two systems, with su ch a kick being ma n- 
dated in the case of XTE J 1118-1-480 (|Fragos et al.ll2007l ). 

An asymmetric kick is not constrained to lie in the or- 
bital plane, as in the case of a Blaauw kick. While the in- 
clination angle of the syst em to the line of si ght is well- 
constrained to be 56° ± 4° (|Shahbaz et al.lll993 ). the longi- 
tude of the ascending node, fi, with respect to the sky plane 
is not known, so we cannot determine the absolute orienta- 
tion of the binary orbital plane. However, for a given value of 
Q., we can determine the orientation of the orbital plane with 
respect to the Galactic axes x, y and z (corresponding to the 
velocity components U , V and W). We assume that the ori- 
entation of the orbital plane does not change with time. For 
any point during its Galactocentric orbit, we can use the 
positions and velocities computed in Section 2] to calculate 
the component of the peculiar motion perpendicular to the 
orbital plane, v± , which provides a lower limit to the asym- 
metric kick velocity should the supernova have occurred at 
that point in time. Assuming all values < f2 < 27r are 
equally likely, we can run Monte Carlo simulations to find 
the probability that the component of the peculiar velocity 
perpendicular to the orbital plane is equal to or less than 
would be expected for the progenitor from the typical ve- 
locity dispersio n of massive st ars (taken as lOkms^^ in one 
component; see lMignardll2000l ) . Fig.[7]shows the probability 
that v± is less than lOkms"^ for the mean and extreme val- 
ues of the distance, as a function of time. The probability is 
low, of order 10-20 per cent, with little variation in the mean 
probability as a function of source distance. It is therefore 
unlikely (although possible, except for certain short periods 
of time) that the measured peculiar velocity perpendicular 
to the binary orbital plane can be attributed to the veloc- 
ity dispersion of the system prior to the supernova. Since 
a symmetric kick cannot give rise to velocity out of the or- 
bital plane, it is probable that there was an asymmetric kick 
during the formation of the black hole. 



6 DISCUSSION 

It appears that a supernova is required to explain the pecu- 
liar velocity of V404 Cyg. Thus the black hole in this system 
did not form via direct collapse. From the range of peculiar 
velocities inferred during the Galactocentric orbit of the sys- 
tem, a symmetric supernova kick could just be sufficient to 
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Figure 7. Probability that the component of the pecuUar veloc- 
ity perpendicular to the orbital plane is less than lOkms^^. We 
assumed a uniform probability for the longitude of the ascend- 
ing node, for angles in the range < fl < 2it, and derived the 
probabilities for the best-fitting distance (4.0 kpc; solid line), and 
the upper (6.0 kpc; dotted line) and lower (2.8 kpc; dashed line) 
limits to the distance. A low probability implies that it is likely 
that the velocity perpendicular to the orbital plane is high, so the 
system has received an asymmetric kick at formation. 



explain the observations if the source is at the lower end of 
the allowed range of distances. However, the mass loss re- 
quired, as well as the component of peculiar velocity inferred 
to lie perpendicular to the orbital plane, make it likely that 
the system was subject to an asymmetric kick during black 
hole formation. 

The full three-dimensional space velocities have been 
measured for a handful of other b lack hole systems. 
XT E J 1118-F480 dMirabel et al.l [iooih and GRO J 1655- 
40 (|Mirabel et al.1 l2002t ) are both in highly eccentric or- 
bits around the Galactic Centre. An asy mmetric na- 
tal k ick is required for XTE J 1118-1-480 (jPragos et al.1 
bOOTl). and believed to have been likely in GRO J 1655-40 
JWiUcms ct~ ai]l2005h. The high-rn ass X-ray binary Cygnus 
X-1 ( Mirabc IfcRodriguesI [2003al ) is moving at only 9 ± 
2kms~^ with respect to its parent association Cyg 0B3, im- 
plying that < IMq was ejected in the natal supernova, and 
that the black hole formed by direct collapse. GRS 1915-1-105 
on the other han d, has been inferred to have an orbit and 
peculiar veloci ty (iDhawan et al.l 20071 ) fairly similar to that 
of V404 Cyg. iDhawan et al.1 l|200'if ). using a Galactocen- 
tric distance of Ro = 8.5 kpc and a circular velocity of 
Oo = 220 km s^^ for the LSR, found that a symmetric su- 
pernova kick and stellar diffusion were insufhcient to explain 
the peculiar velocity unless the source was located at 9- 
10 kpc, where the peculiar velocity is minimized. However, 
using the values of Ro = 8.0 kpc and Go = 236kms~^ as- 
sumed in this paper gives a minimum peculiar velocity of 
23kms~^ for a source distance of 10 kpc, which could then 
in principle be explained by a symmetric supernova explo- 
sion or stellar diffusion. While a distance at the lower end 
of the allowed range would give a higher value for the cur- 
rent peculiar velocity, the peculiar velocity would still be 
sufficiently low at certain points in the Galactocentric orbit 
for an asymmetric kick not to be required. However, if the 
source is at the upper end of the possible distance range 
(> 12 kpc), the peculiar velocity is high enough through- 
out the orbit that an asymmetric kick becomes necessary. 



Without knowing the source distance, we cannot definitively 
determine its formation mechanism. We also note that the 
momentum imparted by a kick of 23 kms~^ to such a 14M0 
black hole (jCreiner et al.|[200lf ) is similar to that gained by 
a ne utron star receiving a kick of a few hundred kms~^ 
(e.g. lLvne fc Lorimeilll994l :lH ansen fc PhinnevI 19971 ) , so ow- 
ing to its large mass, even such a low peculiar velocity for 
GRS 19154-105 would not necessarily rule out a natal kick. 

For V404 Cyg, the derived components of the system ve- 
locity in the Galactic plane, U and V , are much larger than 
W, the velocity out of the plane (Table [21). While this is to 
be expected since the Galactic rotation of 236 km s"'^ forms 
a component of both U and V, but not W, accounting for 
the Galactic rotation (giving the U— < U > and V— <V> 
terms in Table [Sjl does not remove this discrepancy between 
the components of the peculiar velocity in and out of the 
plane. A similar discrepancy is observed for the other four 
black hole systems with measured three-dimensional space 
velocities (iMirabel e t aT '200i|, |2002| : iMirabel fc RodriguesI 
l2003al : loh awan et al.ii2007. ). even accounting for the ranges 
of values allowed by the uncertainties in the system param- 
eters. However, as shown in Fig. O the velocity components 
change with time, so we reconstructed the Galactocentric or- 
bits of all five sources. This showed that the W component of 
velocity can be significantly greater than the peculiar veloc- 
ity in the plane for XTE J 1118-1-480, and can be of a similar 
magnitude to the component in the plane at certain points 
in the orbits of Cygnus X-1 and GRS 1915-^105. With only 
five systems, we are dealing with small number statistics, 
but these results are consistent with there being no preferred 
orientation of the peculiar velocity relative to the Galactic 
plane. Indeed, considering the known population of black 
hole X-ray binaries, including th ose with no measured space 
velocity, the distribution in z (|jonker fc NelemansI l2004h 
demonstrates that a number of systems must have a sig- 
nificant W velocity at certain points in their orbits, in order 
to reach the distances of several hundred parsecs above or 
below the plane at which they are currently observed. Thus 
there is no observational evidence to suggest that the natal 
kick distribution should not be isotropic. 

Thus the two systems with the lowest black hole masses, 
XTE J 1118 +480 and GRO J 1655-40 (both of order 6-7 Mq 
IOroszll2003l ). are found to have required asymmetric kicks 
during the formation of the black hole. For the higher-mass 
systems, the situation is less clear-cut. Cygnus X-1, V404 
Cyg and GRS 1915-^105 all have masses of > IOMq. While 
Cygnus X-1 appears to have formed by direct collapse, the 
peculiar velocity of V404 Cyg implies that a supernova must 
have occurred, and an asymmetric kick seems likely to have 
been required in this system, contrary to asserti ons that 
black holes of lOM© form by direct collapse (e.g. iMirabell 
|2008| ). For GRS 1915-1-105, we cannot definitively determine 
its formation mechanism without an accurate distance to the 
source. Regardless, with only five systems, no clear trends 
with black hole mass can be identified. In order to better 
constrain the formation mechanisms of stellar-mass black 
holes, the space velocities of more such systems with a range 
of black hole masses need to be measured, to constrain the 
frequency of occurrence of natal kicks, and hence supernova 
explosions. 
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7 CONCLUSIONS 

We have measured the proper motion of V404 Cyg using 
20 years' worth of VLA and VLBI radio observations. To- 
gether with the radial velocity and constraints on the dis- 
tance of the system, this translates to a peculiar motion of 
64.1 ± 3. ^lig'skms ^. Given the measured proper motion, 
the black hole cannot have been formed via direct collapse. 
A supernova is required to achieve the observed peculiar ve- 
locity, with either a large amount of mass (~ HMq) being 
lost in the explosion, or, more probably, the system being 
subject to an asymmetric kick. In the case of a pure Blaauw 
kick, ~ IMq must have been transferred from the donor 
to the black hole since the onset of mass transfer, implying 
an initial black hole mass of ~ 9M0 with a donor mass of 
~ 2Mq prior to the onset of mass transfer. 
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